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The stress-energy tensor for the massless spin 1/2 field is numerically computed outside and on 
the event horizons of both charged and uncharged static non-rotating black holes, corresponding 
to the Schwarzschild, Reissner-Nordstrom, and extreme Reissner-Nordstrom solutions of Einstein's 
equations. The field is assumed to be in a thermal state at the black hole temperature. Comparison 
is made between the numerical results and previous analytic approximations for the stress-energy 
tensor in these spacetimes. For the Schwarzschild (charge zero) solution, it is shown that the 
stress-energy differs even in sign from the analytic approximation. For the Reissner-Nordstrom and 
extreme Reissner-Nordstrom solutions, divergences predicted by the analytic approximations are 
shown not to exist. 
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. Classically, the most general static black hole solution to Einstein's equations is spherically symmetric and is 
described by just two parameters, the mass (M) and charge (Q). Quantum effects may be taken into account by 
j3 . solving the semi-classical backreaction equations [1], 

G\xv = 8tt{T^) , (1) 

■ . . 

t-H ' where (T^ v ) is the expectation value of the stress-energy tensor operator, which we shall hereafter refer to as the 
stress-energy tensor for the quantized fields. It includes the effects both of physical particles and vacuum polarization. 
A non-zero stress-energy tensor will alter the black hole's geometry. Even if Eq. (1) is not explicitly solved, finding 
, the stress-energy tensor in a given background spacetime allows one to estimate the size of quantum effects in 
that spacetime. For example, if the stress-energy diverges at the event horizon, quantum effects may prevent the 
existence of the black hole. Such a divergence is predicted in one component of the stress-energy tensor by analytic 
approximations for massless spin 0, 1/2, and 1 fields in Reissner-Nordstrom (0 < \Q\ < M) and extreme Reissner- 
Nordstrom (|Q| = M) spacetimes [2-5] and, in one case, by a two dimensional calculation of the full renormalized 
stress-energy tensor for a massless spin field [6] . 

Numerical computations of the stress-energy tensor for quantized fields in black hole spacetimes provide important 
tests of both analytic approximations and two dimensional calculations. Previous numerical computations of the 
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stress-energy of quantized fields in black hole spacetimes have primarily been for scalar fields in Schwarzschild (Q — 
^Jy 0), Reissner-Nordstrom, and extreme Reissner-Nordstrom spacetimes [7-11,5,12]. While computations relating to 
scalar fields are important and provide some insight into the effects of quantized fields, it is important to also make 
computations with more realistic fields of nonzero spin. Significant differences can exist between fields of zero versus 
nonzero spin [13,14]. Because of difficulties in working with fields of higher spin in curved space, the only previous 
numerical calculation of the full renormalized stress-energy tensor in a black hole spacetime for a field of nonzero spin 
was for the massless spin 1 field in Schwarzschild spacetime [15]. 

In this paper we present the results of numerical computations of the stress-energy tensor for the massless spin 
1/2 field in Schwarzschild, Reissner-Nordstrom, and extreme Reissner-Nordstrom spacetimes in the Hartle-Hawking- 
Israel state [16,17], a thermal state at the black hole temperature. Our results address the validity of the analytic 
approximations that have been derived for these spacetimes [18,2] as well as the question of whether the stress-energy 
is finite on the horizon. We find that in Schwarzschild spacetime the analytic approximation for the field is very poor, 
not even giving the correct sign for the energy density on the horizon. This is in contrast to the cases of massless scalar 
fields [9,10,5] and the spin 1 field in Schwarzschild spacetime [15] where it was found that analytic approximations were 
good both on and everywhere outside of the event horizon. We also find that for Reissner-Nordstrom and extreme 
Reissner-Nordstrom spacetimes all components of the stress-energy tensor are finite on the event horizon. This is in 
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contrast to the analytic approximation of Frolov and Zclnikov [2] which predicts a divergence in one component on 
the horizon. It is also in contrast to the result found by Trivedi [6] that one component of the stress-energy tensor for 
the conformally invariant field in two dimensions is divergent on the horizon of an extreme Reissner-Nordstrom black 
hole. However, it is in agreement with the numerical results found previously for massless scalar fields with arbitrary 
curvature couplings [5,12]. 

In Ref. [19] a method was given to compute the stress-energy tensor for the massless spin 1/2 field in a general 
static spherically symmetric spacetime. The method involves writing the stress-energy tensor in terms of the Euclidean 
Green's function for the field and renormalizing using the method of point splitting [20,21]. An explicit expression for 
the Euclidean Green's function in a static spherically symmetric spacetime was derived in terms of mode solutions to 
the Euclidean Dirac equation. The time dependent part of the modes and the angular dependent part were both given 
in terms of known functions. The radial part was given in terms of solutions to a coupled set of ordinary differential 
equations. For most spacetimes, including black hole spacetimes, these equations must be solved numerically. The 
final expression for the stress-energy tensor consists of two parts both of which are finite and separately conserved. 
For one part the mode sums have been computed analytically for a general static spherically symmetric spacetime. 
It has a trace equal to the trace anomaly and thus can serve as an analytic approximation for the entire stress-energy 
tensor. It is equivalent to previous analytic approximations for the massless spin 1/2 field [18,2]. The other part must 
usually be computed numerically. 

To simplify the numerical computation of the stress-energy tensor it is helpful to use the WKB approximation for 
the high frequency and/or large angular momentum modes [9,10,5]. The higher the order of the WKB approximation 
used the fewer modes that must be numerically computed to give an answer to a given accuracy. However the number 
of terms in the WKB approximation grows rapidly with order so it is not practical to use too high an order. 

The metric for a Reissner-Nordstrom spacetime can be written 

ds 2 = -f(r)dt 2 + f{r)- l dr 2 + rW , (2) 

where 



The temperature of the black hole is 



2nT = (4) 

2M 2 -Q 2 + 2M V / M 2 - Q 2 

Recall that in the Hartle-Hawking-Israel state, the field is assumed to be in a thermal state at this temperature. Note 
that the extreme Reissner-Nordstrom black hole is at zero temperature. 

In Figs. 1 and 2 our results for Schwarzschild spacetime are displayed. As in the case of conformally invariant spin 
[9,10] and 1 [15] fields, the weak energy condition is violated near the horizon 1 . The analytic approximation is very 
poor near the horizon (although it is accurate far from it), and even gets the sign of the energy density wrong. This 
is in contrast to the spin and 1 cases where the analytic approximation is good near the horizon [9,10,15]. 

As a check on this surprising result, we also computed the energy density at the horizon using radial point splitting 
with one point fixed on the horizon. This method was developed by Candelas [7] who used it to compute the 
stress-energy tensor at the horizon for the conformally invariant scalar field in Schwarzschild spacetime. It has the 
advantage that modes of only the two smallest allowed frequencies contribute to the mode sum. The method in [19] 
involves splitting the points in the time but not the space direction. The resulting expression for the stress-energy 
tensor contains a sum over modes with arbitrarily large frequencies. Thus the two methods, while not completely 
independent, arc significantly different. The renormalization counterterms are also different for these two methods. 
The energy density of the spin 1/2 field on the horizon in the radial point splitting calculation was found to agree to 
better than four digits with that obtained in the time point splitting calculation. 

Some of our results for Reissner-Nordstrom spacetimes are displayed in Fig. 3. For non-zero charge Q, there is a 
well-known ambiguity in the analytic portion of the stress-energy tensor caused by renormalization [21,5,19]. This 
ambiguity docs not arise for the Schwarzschild metric, nor does it exist on the horizon for (T* t ), (T r r ), or (T e g) in 



1 With our conventions the energy density is p — —T t *. 
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Reissner-Nordstrom or extreme Reissner-Nordstrom spacetimes. For definiteness, we have chosen the renormalization 
parameter /i so that fiM = 1. From Fig. 3, it is clear that the energy density on the horizon is negative and its 
magnitude increases with increasing \Q\/M. Thus the weak energy condition is violated in these cases as well. 

One of the predictions of the analytic approximation for the massless spin 1/2 field is that one component of the 
stress-energy diverges on the event horizon for all Reissner-Nordstrom black holes. The component is related to the 
energy density seen by a freely falling observer passing through the event horizon and is explicitly 

W) • (5) 

For \Q\ < M the divergence in this component of the analytic approximation is logarithmic and of the form 

Q 



207r 2 r 6 



log[ M 2 M 2 /M] , (6) 



where /x is the renormalization parameter discussed above. This divergence also occurs in the case \Q\ = M . However 
in that case there is a term in the analytic approximation for the component (5) with a stronger divergence. It is 

M 2 (32M-23r) 

1807r 2 r 6 (r-M) ' ( ' 

This is also the type of divergence found by Trivedi [6] for the same component of the full renormalizcd stress-energy 
tensor of a conformally invariant scalar field in a two dimensional extreme Reissner-Nordstrom spacctime. 

However, when we numerically computed the component (5) for the massless spin 1/2 field in Reissner-Nordstrom 
and extreme Reissner-Nordstrom spacetimes, we found no evidence that this component of the stress-energy tensor 
diverges on the event horizon. This is in agreement with results found for massless scalar fields with arbitrary curvature 
couplings in these spacetimes [5,12]. Some of our results are shown in Fig. 4. 

In summary the full renormalized stress-energy tensor for the massless spin 1/2 field has been numerically computed 
in Schwarzschild, Reissner-Nordstrom, and extreme Reissner-Nordstrom spacetimes in the case that the field is in the 
Hartle-Hawking-Israel state. The analytic approximation for the stress-energy tensor in Schwarzschild spacetime is 
found to be very poor near the event horizon. Divergences in the stress-energy, predicted by the analytic approxi- 
mations for Reissner-Nordstrom and extreme Reissner-Nordstrom, are found not to occur. Thus it is clear that the 
analytic approximations for the massless spin 1/2 field in static black hole spacetimes are of extremely limited use 
and must usually be supplemented by numerical calculations of the full renormalizcd stress-energy tensor. 
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Figure Captions 

FIG. 1. The analytic (dashed) and numerical (solid) values for the quantity (T* t ) as a function of distance from the event 
horizon r+. 



FIG. 2. The numerical values for the quantities (T 6 ' e) and (T r r ) from top to bottom at the event horizon r+. 



FIG. 3. The numerical values for the quantity (T* t ) for \Q\/M = 1, 0.99, 0.95, 0.8, from top to bottom at the event 
horizon r+ for the case jiM = 1. 



FIG. 4. The numerical values for the quantity ((T r r ) - {T\))/f for \Q\/M = 1, 0.9999, 0.999, 0.99 from top to bottom at 
the event horizon r+. For these plots, the value \xM = 1 has been chosen. 
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